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APPLICABILITY OF THE AVAILABLE SCSiPTION 
MODELS FOR PREDICTING MERCURY REMOVAL 

ABSTRACT 

All the investigations r^orted todate in the area o£ heavy 
metal removal from water were concerned mainly with giross removal 
efficiencies and no attempts were nade to formulate any generalized 
mathematical approach for predicting the performance of sorption 
systems. On the contrary, the organics-sorption understanding has 
progressed considerably during the last decade and many mathematical 
models have been proposed by different workers to predict perform- 
ance characteristics of sorption systems. In the present study, a 
caipr^ensive listing of the available sorption models was prepared 
and on the basis of the available data two modified approaches were 
proposed to predict the effluent concentration- time profiles for 
mercury-coal and mercury -active carbon systems. The Ihomas solutim 
for the second-order kinetics supplemented by Hiester ar^ Verm^len, 
using theoretical as well as graphical methods for the evaluaticai of 
the sorption rate constant, constitutal the first approach. In the 
secoi*"' approach the original lumped parameter method used by 
Keinath was irodified into a nondimensional form and the equaticais 
weire solved numerically by the fourth order Runga-Kutta method. 
Weighing the modified approaches with the limited data available 
revealed that the first approach could reasonably predict the per- 
formance of mercury sorption systems with influent oohcentration in 
the milligram range, The second approach, i.e., the modified lumped 
parameter approach predicted a pretty early breakthrcaagh. The limi- 
tations of the two proposed approach^ were discaissed in light of 
the difficulties involved in reli^le estimation of the dSi.fereit 
model parameters. 
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1, INTRODUCTION 


After the Minamata mercury pollution episode in Japan 
in 1953, heavy metal pollution began to concern the guardians 
of public health. Their possible presence in the ecosystem 
and the environmental effects of mercury, lead, cadmium, 
chronixom, silver, etc, were investigated. The analytical 
methods including instrumental analysis for detecting minute 
concentrations of these metals were studied and improved upon. 
Analysis of low concentrations of heavy metal compounds is 
very difficult but significant even for fixing up their 
permissible limits in drinking water. Atomic absorption 
spectrophotonetry is recognised presently as the most reliable 
of all the proposed methods (Standard Methods, 1975), Once 
the analysis part was performed confidently the environmental 
engineers started investigating the efficacies of the conven- 
tional water and wastewater treatment processes in removing 
ho vV metals (Logsdon and Symons, 1973), 

Tue first major challenge to the environmental 
engineers was thrown by the epidemic diseases which was pretty 
well tacdcled by including the chlorination step with the 
traditional treatment processes (McCarty, 1980). The second 
challenge x-ta-s thrown by the organics, This time too they 
emerged winners by virtue of the grand succsss of active carboi 
sorpticn of organics. If heavy metal pollution is considered 
to be the next big challenge, it should be admitted that the 
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fight is still continuing as workers in this field are busy 
studying the removal efficiencies and mechanians of removal by 
the different processes. 

It has been observed that active carbon whidi played 
the tricSc for organics r®noval, does also adsorb heavy metals 
hut the extent of removal is definitely lesser. In the case 
of mercury ronoval, this led seme workers to opt for cheaper 
sorbents like wool, hair, rubber, coal etc, (Michelsen and 
co-workers, 1975; Pandey and Chaudhuri, 1981), But all these 
studies on heavy metal removal by sorption are rather quali- 
tative and limited to specific systems only. On the other 
hand, the sophistication of organics -sorpticai research has 
increased by leaps and bounds during the last decade and a 
number of complex mathematical approaches have been proposed 
(Crittenden and Weber, 1978a; Keinath and Weber, 1968; Keinath, 
1975; McKay, ottertaum and Sweeney, 1980; Wilson and Clarke, 
1979), most of them being solved numerically with modem 
computers. 

To keep pace with the progress in organics -sorption 
research it was felt that the available sorption models should 
be oempixed and their potential applicability for predicting 
heavy metal renoval should be investigated. It would be 
logical, as a first step, to include the pure, single-solute 
models to be follov?ed by the mom complex models developed for 
bi- and multi -solute rystems. Such an approach should ulti- 
mately lead to a more appropriate model taking into consider- 
ation the limitations of the existing models developed for 
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organics sorption. The present study was limited to the first 
step with particular reference to mercury. Ihie models are 
listed here and the ones developed for pure^ single-solute 
systems are assessed for application to the mercury sorption 
system ultimately leading to seme modified approadhes for 
final weighing using the data reported by Humenick and 
Sdnnocr ( 1974) and Pandey ( 1978) for mercury- active carbon 
and mercury-coal systems, respectively. 
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2, AVAILABLE SORPTION MODELS 

During the last decade numerous models have been 
proposed to explain sorption o£ organics and dyes onto active 
carbon. These models range from mathematical and conceptual 
descriptions of breakthrough curve and numeric solution models 
(Crittenden and Weber, 1978a) to detailed investigations of 
the thermodynamic properties of the solvents, sorbents and 
sorbates (McGuire and Suffet, 1978). The models are not 
equivalent in data needs, theoretical justifications or appli- 
cations, but are useful as a basis for understanding the 
sorption mechanism, A complete list of the available models 
is presented here. It is interesting to note, however, that 
till now none of these models has been applied to heavy metal 
sorption. The discussion following the list also includes 
preliminary assessment of the models to qualify their selec- 
tion for a more detailed weighing using the avail^le data 
on mercury sorption. 

Model Reference 

A. Old Models Swan and Urquhart (1927) 

a. Lagergren Rate Equation (1898) 

b, Arrhenius Equaticai (1910) 

c. Pickles Equation (1920) 

d, Gustaver Equation (1922) 
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Model Reference 

B. Conceptual Models ' McGuire and Suffet (1973) 

a. Mass Transf r Zone 

b. Ideal Adsorbed Solution (IT. 3) 

Theory 

c. Net Adsorption Energy Concept 

C. Michigan Adsorption Design and Crittenden and Weber 

Appl icaticais Model (MADAM)' (1978a) 

D. Logistic Curve oulman (1980) 

E. Diffusion Model McKay, otterburn and 

Sweeney (1980) 

Keinath (1975) 

Wilson and Clarke (1979) 

Thomas (1944) 

Miyamoto and Tram (1979) 

Famularo, Mueller and 
Pennu (1980) 

K. Convection Limited Mcjdel Hendricks (1980) 

A, Old Models 

Swan and Urquhart (1927) rave e very comprehensive 
review of the cxjntemporary sorption models available in the 
first q.uarter of the present century. Out of a very long 
list given, the following four models v7ore subjeerted to a 
preliminary assessment using tlae available data cn mercxiry 
removal reported by Eumenick and Schnoor (1974; and Pandey 
(1978). 

a, uagergren Rato Eejuation 

= k(x^^ “ X. ) (i) 

eejra t: 


F. Keinath Model 

G. Wilson and Clarke Model 

H. Second Order Kinetics Using 
Themas Solution 

I. Miyamoto and Trem Model 

J. Famularo Concept 


dx 

5t 
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where, k = rate cf sorption, 

X = mass of sorbate sorbed by unit weight 
eqm 

of sorbent at equilibrium', and 
= mass of sorbate sorbed by unit wei^t 
of sorbent at time, t. 

This equation does not conform to general requiroments of the 
rate equation and Gustaver (Swan and Urquhart, 1927) had pointed 
out that it is inapplicable to high concentrations, since it 
assumes constancy of concentrations during sorption. 

b, Arrhenius Equaticn 

log - 0.4343 | = kC (2) 

S ^ X 

where, C = equilibrium ccxicxaitration of solution, 

s = mass of sorbate sorbed per unit weight of 
sorbent at saturation, and 
X = mass of sorbate sorbed per unit weight of 
sorbent. 

c. Pickles Equation 


k 


1 

^2 ^ H 


a 

log - 



0,4343 



(3) 


where, a = original amount of sorbate, aid 

^l'^2 ~ ^^^ss of sorbate sorbed per unit wei^t 
of sorbent at tiit^ t^, t: 2 * 


This ©juaticai had been used for -sorption from gas as well as 
fran liquid, 

d. Gustaver Equaticxi 



tla 


1 

■ X ) 

eqm 


log( 


a .. 

5 ^ 

eqm 


-|3S1) 

a 


(4) 


X 
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Gustaver obtained the equation hy assuming that the amount of 
sorptican was proporticsial to the free surface and the total 
amount of sorbate in ihe liquid jhase* 

All the above four models are empirical and applicable 
to specific systems only. They were found unsuitable for 
application to the mercury system on the basis of preliminary 
assessment. Ihe detailed analysis, however, is not presented 
here. 

B. Conceptual Models 

a. Mass Transfer Zone 

The height and rate of movement of the mass transfer 
zone in an active carbon column was used to characterize the 
effects of the process variables on sorption of organics 
(McGuire and Suffet, 1978). The height of the mass tr^sfer 
zone is a measure of the rate of removal of the solute by 
the sorbent and the rate of movemoit of the mass transfer 
zone is a measure of the affinity of the solute for the sorbent 
surface. The model is applicable when the solute has a rela- 
tively high affinity for the sorbent surface. Attempts to 
use the mass transfer zone concept to define the breakthrough 
characteristics of poorly sorbed organics were unsuccessful. 
As heavy metals do not have a very high affinity for active 
carbon this model was not considered any further in this 
study. 

b. Ideal Adsorbed Solution (lAs) Theory 

This theory describes a correlation for estimating 
the sorption of organic compounds faxm dilute aqueous solutions 
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onto active carbon (McGuire and Suffet, 1978). A formal 
statistical.<tiechanical treatment of sorption from dilute 
solutions provided the basic equation to develop the corre- 
lation, Developing the correlation for mercury-sorbent systems 
was considered to be beyond the scope of the present study, 
c. Net Adsorption Energy Concept 

The calculated net adsorption energy is a function 
of a compound’s aqueous solubility, molecular weight, density, 
dispersion cmponent of the solubility parameter for the 
organic cxmpound and that for the active carbon surface. 

This concept was used to evaluate sorption data for medium- 
polarity to high-polarity organic compounds and to relate that 
to other sorption data in the literature so that a general 
understanding of sorption could be approached through an 
understanding of cxxnpounds and surface polarities (McGuire and 
Suffet, 1978), For the mercury -active carbon or mercury-coal 
systems this understanding ^out the compound and surface 
polarities are not very clear and hence the net adsorption 
en^^i.gy concept cannot be applied at present, 

C. Michigan Adsorption Design and Applications Model (MADAM) 
Crittenden and Weber (1978a) developed the Michigan 
Design and Applications Model (mAdAM) to facilitate perfor- 
mance. forecasting for optimum design of sorption systans and 
for use as a predictive ccanponent of a feed-forward operational 
contix)! system. They described conceptual and mathematical 
bases underlying the model , techniques used for parameter 
estimation and methodologies for determination of surface 
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diffusion coefficients and film transfer coefficients* 'Iliey 
also tested the models for single- and multi-component 
solutions using independently determined model parameters 
(Crittenden and Weber, 1978b), It was projected that complex 
factors sudi as fluid dispersion, solidsHtxixing, multi-solute 
Interactions and surface biological growths, could eventually 
be incorporated into the numeric solutions. The partial 
differential equation for licya id-phase material balance was 
simplified and the numeric method, the implicit finite back- 
ward difference method was compared by Weber and Crittenden 
(1975) to the famous analytical solution proposed by Thomas 
(1944), They found that the solutions gave almost identical 
traces when plotted and that the error was much less than 
one percent. For solution of solid-phase partial differential 
equation, the Crank-Nicholson method was selected, Ihe influent 
and effluent concentration-tirr.s profiles were given for phenol, 
p-brcmophenol, p- toluene sulfonate and dodecyl benzene 
sulfonate (Crittenden and Weber, 1978b). 

The solutions provided in the publicatiois referred to 
above, are not non-dimensional. So these solutions could 
not be applied for mercury-active carbon or mercury-coal 
system. The computer programs for the solution have not been 
published. So if one wants to apply the MAn^M for mercury 
sorption then one would have to solve the liquid- and solid- 
phase partial differential equations by the implicit finite 
backward differencse and Crank-Nicholson methods respectively, 
which was consideraS to be beyond the scope of the present work. 
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D. Logistic Qirve 

The logistic curve was proposed as a model for deter- 
mining the design parameters for active carbon coltunns 
(Oulman/ 1980), The parameters were used for scaling up the 
design of a carbon column fran laboratory and pilot scale to 
full scale, 

Ihe Bohart-Adams equation was published in 1920 for 
predicting the service life of the active carbon used in a 
gas mask for the removal of chlorine fran air. The equation 
can be written as 

InizS' - 1) = K'NxVv - K‘C t (5) 

u a o 

where/ C = effluent concentration/ 

C = influent concentration/ 
o 

K‘ s= sorption rate coefficient, 

N = sorption capacity coefficient, 

= bed depth/ 

v = approach velocity (rate of applicaticai), and 
t = service time. 

Equation (5) can be rewritten as 



(a + bt) 


1 + e 


( 6 ) 


The plot of C/C^ vs, t is the Logistic Curve which is S -shaped. 
The data from a sorption run can also be linearized 
according to equation (5) by plotting ln{C/c^ - 1) vs, t, 
Althtxigh this is a perfectly satisfactory way of handling the 
data, biologists have been using the inverse expression for 
many years in linearizing dose response relationship 
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ln( 


C/C 


1 - C/C, 


- K'Nx,A^ + K’C t 
a o 


(7) 


The left hand side of equation (7) is called the "logit” of 
C/C^ whidi is the contraction of the term "logistical unit". 
Frcxn the plot of logit of C/C^ against t, average values of K 
and N can be calculated which may be used in predicting full 
scale plant performance. 

The column studies on mercury sorption reported ty 
Humenick and Schnoor (1974) and Pandey (1978) were not conti- 
nued upto breakthrough. So logit of C/c^ remained almost a 
constant and it was not possible to calculate K and N, the 
sorption rate and capacity coefficients, 

E, Diffusion Model 

A s« 2 hool of thought (Mcacay, otterbum and Sweeney, 1980) 
believes that the mechanian of sorption involves a study of 
the diffusional processes. Pick’s second law of diffusion 
in terms of spherical polar co-ordinates r, © and 0 is given 


by (Crank, 1956) 


dC 

3t 


T Tr + ^ ^ ^ 


1 

"in© r'l© 


a. 

■T© + 




s?n © o0 




where, D = molecular diffusivity. 

The following assumptions were made for solving the equaticai 
(Skelland, 1974) : (i) The crsnoentration of solute is uniform 
at C^ throughout the solution and is zero in the sorbent 
particle at the start of diffusion (t = O), (ii) The resistance 
to transfer in the mediiam surrounding is negligible, and 
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(iii) Diffusion is radial/ there being no variation in caancen— 
tration with angular positiai. The equation was solved using 
the initial assumption which specify the boundary conditions 
as C = C and C = where, r = particle radius, ^d 
a= equilibrium concentration. The fraction solute sorption 
onto the sphere at time t is given in equation (9) where C 
is the average solute concentration throughout the sjhere at 
time t 



1 - 


6 

3 


00 




1 

n 


exp 


2 2 

/Dn'^ ic^tN 

( S5 — } 


(9) 


where, n = an integer, 

A numerical solution is found in the literature (Newman, 

1931) in the form of a plot of 1 - (C - C)/(c - ) vs, 

o ' o oo 

2 

Dt/r (Pig, 1), The solution is non-dimensional and is, 
therefore, apparently very useful for application to the 
mercury systems under consideration. But the crux of the 
problem, which causes the limitation to this model, lies in 
the difficulty in determining the actual value of the diffu- 
sivity. When intraparticle diffusion is accomp^ied by 
sorption, as in the case of the mercury systems, evaluation of 
characttristic diffusivities from experimental rate data can 
be accomplished cxily by appropriate modification of Pick's law 
for diffusion to include a term for the sorption process (Weber, 
1972). If the relationship between the concentration of 
sorbed solute and sc lute free to diffuse is linear, the 
modification of Pick's law still permits rather steady analytic 
solution for the characteristic diff us ivity. If the 
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FIG. 1 UNSTEADY-STATE MOLECULAR DIFFUSION 
IN A SPHERE, SLAB OR CYLINDER 
(Newman, 1931). 







i 

1 
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relationship between free and sorbed solute is nonlinear, as 
in the mercury systems, the resulting modification of the 
expression for Pick's lav.* is nonlinear and numerical integr- 
ation methods must be used to obtain a solution for the 
characteristic diffusion coefficient (Weber, 1972). 

F, Keinath Model 

Keinath (1975) published a quite detailed discussion of 
the modeling of sorption contactors by means of a lumped 
parameter approach and included computer programs for carrying 
out these calculations. 

From material balance relationship for a packed-bed 
case, the developed partial differential equation upon simpli- 
fication becomes 

TT is 4 . ^ ^ As _ o ' (10) 

where, U = solution volur-''tric flow rate, 1/hr, 

C = solution-phase concentration of solute, mole/1, 
q s= solid-phase concentration of solute, mole/g, 
z = cKial distance, an, 
e = void ratio, and, , 

t = time, hr. 

Keinath (1975), instead of solving this equation, adapted the 
lumpsed parameter approach through the use of a continuous 
system modeling program, CSMP/360, produced by 135. Solutions 
are given for four cases, viz,, i) single— solute syston, 
constant input, ii) single-solute syston, sine-curve flow 
input, iii) single-solute, countercurrent system, and 
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iv) di-solute system. As the available data for mercury- 
sorbent systems are for the single-solute system with ccnstait 
input/ only the first case is discussed here. 

Use of the lumped parameter approach required segmen- 
tation of the packed-bed of sorbent into a discrete number of 
finite elements {Friedly, 1972). Keinath divided his column 
length (about 32,4 cm) into ten such elements. Both the 
liquid- and solid-jhase concoitrations of solute were assumed 
to be uniform throughout each element. Maintaining ccaitinuity 
of mass flow of the solute between adjoining elements was 
also important for both the liquid and solid phases, Keinath 
(1975) proposed the following formulation for an element of 
the packed-bed 


Solution -phase 


Input to element 
Output frcxn element 
Sorptiai 
Accumulaticai 


^^n-1 

^^n 

e (dC^/dt)v 


Solid- ihase 
O 
O 

P^A^ 

f (dq^/dt )V 


Where, V = volume of elements, liters, 

= rate of sorption, raole/g-hr, and 
n = number of elements. 

From material balaice: 

Input to element = Output f ran element + Sorption 

+ Accumuiaticm 

The final equations 


dC TT 

a# = «iv><<=n-l - ^n^ 


(^/e )R, 


( 11 ) 
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& - 

itiuS/ adapting the lumped parameter approach facilitated 
elimination of the partial differential equation. Ihe 
essence of Keinath approach was th'- development of equations 
(11) ^d (12) and solving them ;n a straightforward manner by 
using CSMP/360, 

If the CSMP/360 program w’ero available, application of 
Keinath model for the mercury systems would have been easier. 
Since this program was not available in the dec- 10 system at 
IIT-Kanpur, direct application of Keinath approacdi was not 
possible. In spite of this setbadc, attempts were made to 
solve the equations by the fourth order Runga-Kutta method 
(Chapter 4), very standard numerical method for solving 
ordinary differential equations, !7ith a view to making a 
general approach, the equations xvero modified into their non- 
dimensional form (C.apter 4), 

G. Wilson and Clarke Model 

This model v/as also based on cho lumped parameter 
approach like the Keinath model. Sorption of solutes onto 
active carbon was modeled at the microscopic levels by means 
of the equations of continuity and mass balance. Two competing 
solutes were assumed to be present. The results of this 
approach were then cempared with the much simpler lumped 
parameter model for active carbon sorpticai for two competing 
solutes. Quite good agreement was obtained validating the 
utility of the simpler model (Wilson and Clarke, 1979; 

Wilson, 1979). 
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Tfie available data for the mercury systems being for 
pure, single-solute systems only, the Wilson and Clarke model 
framed for bi- solute systems was not included in the present 
work. As an example of the lumped parameter approach, only 
the Keinath model was studied in details, 

H. Second-order Kinetics Using Thomas Solution 

The rate equation (13) and the conservation of mass 
equaticffi (14) together provide a description of the mass- 
transfer process sufficient to permit a mathematical solution 
(Weber, 1972) 


d(q/q^) 


dT 


C_ 

C 


(1 


^) 

^oo 


r (1 

‘^co 


•) 


t 

s 


C (V 
o 


vfg ) 


‘oo 


eb ^ 


(13) 


(14) 


where, r - equilibrium parameter = ■ ■ ■ ^ 


s = column capacity parameter = 


k g KSL 


= solution capacity parameter = 
kf. 


'E 


1 ^ 
% + 


) 


V = bulk volume of column, under consideration, 
v^ = macroscopic water flow velocity, 

V = total volume of feed that has entered the 

column , 

= det sity of dry sorbent in column, 

f„ = fractitXi of external voids per unit volume 

■jE. ■ , 

of the sorbent packed. 
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k = sorpticai rate constant, 

K = Langrraiir isotherm constant, 

L = length of column, 

R = volimietric flow rate of fluid through fixed 
sorbent, and 

S = ultimate capacity of dry sorbent. 

These equations were integrated, subject to appropriate, 
boundary conditions, by a method described by Thomas (1944) to 
yield the following expression for the ratio of the concen- 
traticais of solute in the effluent stream, C, to the initial 
concentration in the influent stream, C^, for any specified 
combination of the parameters r, s and t 

C J(rs,t) nci 

J{rs,t) + [i - j{s,rt)] exp [(r-l)(t-s) ] 

X 

whei?e, J(x,y) = 1 - e”^ 5 y''^)dz, 

o 

X, y and z *= integral variolas, and 

I, s= modified Bessel function of tiie 

s- 

first kind. 

Among the usable mathonatical approaches for prediction 
of efiZucnt concentration in a column syst®n, the method of 
Thonas (1944) has the greatest generality, because only this 
method takes into account the full range of possible variations 
in equilibrium constats for sorptiixi reactions (Hiester and 
Vermeulen, 1952; Vermeulen, 1958; Vermeulen, Klein and Hiester, 
1973; Vilker and co-workers^ 1978). The theory of Thonas 
has escaped the use it deserves mainly because of the functions 
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which are apparently difficult to evaluate, Hiester and 
Vermeulen (1952) gave the following workable form of the 
expression originally developed by Thanas (1944) 



^ f l - H(J^ - /s)l exp( ~ r l/( 4- ^'^■/rst)l 

/i [l - H(vrs “ '/t)] exp( vrs - '/t) + [l/(»/t + ^Jrst),f] ’ 

(16) 

where, H(z) = error function. 

The error functions may be noted fran the tabulated value and 
the value of the whole expression may be calculated using a 
pocket calculator. Hiester and Venneulen (1952) also published 
plots of C/Cq against t/s for r = 0,0, 0,5, 1.0, 2,0, 10,0 
and 100,0 (Graph for r *= 1 is shown in Fig, 2), 

In fine, the procedure for predicting the concentration 
profile for a column system is as follows 

i) r, s and t are calculated from their known expressions, 

ii) if r has any of the values for which graphs are available, 
the appropriate graph is taken and the C/C^ values are 
predicted for the known values of s and t/s (Pig, 2), and 

iii) if r has some other value, the expression (16) is 

directly evaluated by using the tabulated values of error 
functions. 

Looking carefully into the ^ove procedure it is obvious 
that as far as calculations are concerned, the crux of the 
problem lies in finding true values of the parameters s and t. 

If one looks back into the expressions for s and t 
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THROUGHPUT PARAMETER (f/s) 


FIG. 2 GENERALISED BREAKTR0U6H CURVES 
(Hiastar and Varmaulan, 19S2) . 
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(17) 

(18) 


one finds that all the terms except k are either known or may 
be found out experimentally. The available expression for k 
is as follows (Hiester and Vermeulen^ 1952) 


1 

Ic 


1 + (r“l)(q/q ) (l-r)(C/C ) + r 

oo ^ o 






(19) 


where, kp = mass transfer coefficient for case of fluid- 
liiase transfer controlling, 
kp = mass transfer coefficient for case of transfer 
inside solid particles, 

ap,ap = effective areas for mass transfer for tiie 
respective cases, and 

q ^ 

Dq = ;r — 5 — = ratio of concentrations on solid- 
^o E 

and in liquid-^Aiases at saturaticxi. 

So now for finding the true values of the parameters s ^nd t, 
one would have to determine the transfer coefficients and 
the effective areas. The determination of the transfer coeff- 
icients, even though very difficult, is facilitated 1 ^ the 
available expressions (Keinath, 1975) 

60 D 

= — - P - (1 - fE> 

Where, Op = particle diameter. 
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fg = porosity, and 

D^re - pore diffusivity term. 

Again, determination of the value of the pore diffusivity term 
is difficult 



where, t= diffusivity of solute in solvent, and 
fj = internal porosity. 

The other mass transfer coefficient kp may be determined frtxn 
expressions containing Reynolds nimber, Sherwood number and Schmidt 
number. This method of finding kp was used by Keinath (1975). 

The other method (Famularo, Mueller and Pannu, 1980) of 


finding kp is based on the following equation 



V 


( 22 ) 


where, V = batch volume of liquid. 

At the beginning of a batch rate run when the solid-phase 
concentration is very small, it can be shown that liquid 
resistance clearly controls the overall transport of solute 
into the particle (Famularo, Mueller and Pannu, 1980), Under 
this cx»ndition only the equation (22) is applicable for deter- 
mining kp. 

Coming back to the problem of evaluating the parameters 
s and t, one may comment, perhaps without being pragmatic, 
that these very sensitive parameters can never be determined 
without reservations and hence the potential of Thcmas solution 
is reduced significantly by this weakest link. 



23 


I. Miyamoto and Tram Model 

Ito eliminate the just discussed indirect method of 
determining the sorption rate c»nstant, k, Miyamoto and Tram 
(1979) proposed a direct graphical method. They studied the 
sorption of a dispersed dye by sand and calcareous soil using 
laboratory batch and cxilumn bests for the purpose of obtaining 
cjuantitative information needed for disposing of dyeing 
effluent. Miyamoto and Tram (1979) proposed tlie following 
rate equation for sorption data fitting well with the Laignuair 
isotherm 

-~ = k [K(S - s)C » s] (23) 


where / k = sorption rate constant/ 

. S rs ultimate capacity of dry adsorbent, 

K = Langmuir isotherm ccaistant, 
s = amount of sorbate sorbed, and 
C = equilibrium dye cxaicentration. 

At equilibrium, =0. Thus, the right hand side of 
ecjuation (23) reduces to 

s « K(S - s)C 


(24) 


The eciuotion (24) is nothing twt the representation of Langmuir 
sorption isotherm. Fran the law of mass conservation 

ic 

a t 


- ^ s 
= - ® Tt 


es(t) = - c(t) 


(25) 

(26) 



24 


where ^ e =: soil/solution ratio. 

Using equations (25) and (26), the equation (23) was rewritten 
as 

= - k [kc^ + (e KS - KC^ + 1)C - C^] (27) 

Integrating equation (27) from t = O to t = t and frcm C = C 

o 

to c = c 

ln( 2KC^ + a + p ) ( 2KC + a - p )/( 2KC^ + a - p ) ( 2KC + a + p ) 

= “p kt (28) 

or In X = - P kt 

where, a - e KS - KC^ +1, and 
p .= (4KC^ +a^)^/2. 

If the plot of (- 1/p )ln X against t is linear (Fig. 4a) then 
the rate constant k can be determined fron the slope of the 
plot* The rate constant was determined for mercury-coal 
systan and was used in Thcmas solution to determine the 
column and solution capacity parameters s c«id t (Chapter 4). 

J. Famularo Concept 

Fatnularo, Mueller and Pannu (1980) conceptually descri- 
bed the oaibon models, viz,, homogeneous solid diffusion 
model, pore diffusion model and uniform state model. They 
proposed a micropore-raacroshell particle concept. The decline 
in concentration, according to them could roughly be divided 
into an initial high- rate interval followed by a lower rate 
approac^a to equilibrium. Such behavior is consistent with the 
porous nature of active carbon and can be attributed to a high 
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rate of solute movement through raacropores followed by mucii 
slower passage tharough micropores. The detailed model 
(Famularo# Mueller and Pannu^ 1980) was not consideared in this 
study as it involves several mass transfer coefficients and 
other parameters which are difficult to determine for mercury- 
sotbent systems, 

K. Convection Limited Model 

The ancillary kinetic models focus upon the various 
diffusion mechanisms at the particle level, e.g.# film diffusion, 
pore diffusion and pore-surface diffusion. Another sorption 
rate mechanism not considered thus far is the bulk transport of 
sorbate to the particle by convection (Hendricks, 1980), The 
convection limited model supplements the established theory 
based upon the material balance equation and particle kinetics. 
The theory is applicable particularly to the design of packed- 
bed sorption reactors in which the sorbate feed concentrations 
are very low and this model should be applicable for removal 
of low concentrations of heavy metals from drinking water. 

For low concentrations of sorbate, if the sorption rate is 
calculated on the basis of diffusion only the predicted rate 
would be higher than the actual rate, 

Iho available mercury-sorption data were not suff icLait 
for testing this model and hence the potential of this very 
interesting mcjdel remained undetermined at this stage. 



26 


3. SPECIAL AS>ECTS RELATED TO MERCURY SC21PTION 


This chapter includes a discussion on scxne special 
aspects that need to be included or taken into consideration 
while attempting to predict mercury removal through the use of 
mathematical modeling. A section on the possible hazards of 
active carbon or charcoal as sorbent in water treatment is 
also included . 

A, Species Distribution of Mercury Compounds in Water 

When Hg and Cl are present in water they will combine 
to form the undissociated but soluble species of HgCl 2 . 
Chloride can also coabine with mercury in other proportions 
to form a variety of complexes and equilibrium relationships 
can be developed for the various mercury -chloride species 
from the following reactions (Sawyer and McCarty, 1978) 

= 5.6 X 10^ 

K 2 = 3.0 X 10 ^ 

K 3 = 7.1 

K4 = 10,0 


Hg++ + cr 


-Hgcl 

HgCl"^ + Cl" 


H^l, 

HgCl2 + Cl" 


HgCl' 

HgCl3 + Cl" 


HgCl 


Reimers, Krenkel and Eagle (1974) gave a table showing the 
species distribution of tiie inorganic foirms of mercury under 
different conditions of chloride concentrations and jSH, 
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Chloride Concentration Species of Inorganic Meroiry 



s= 5 

pH » 7 

pH = 9 

18000 ppm 

HgCl** 

HgClJ" 

HgO,H20 

10000 ppm 

HgCl““ 

HgCl’" 

Hg0.H20 

100 ppm 

Hgci2 

HgO,H20 

Hg0,H20 


According to Morris and Weber ( 1966) , intraparticle 
transport phenomena may derive from the net action of several 
molecular forces. For most sy stasis of interest in water and 
waste treatment applications, these forces may be enianerated 
as sorption, two and three dimensional micelle formation, 
electrokinetic interactions and molecular diffusion (Weber, 
1972). In the studies of organics-sorption, the electro- 
kinetic interacticn was neglected in general. But in the case 
of mercury -sorption, from the preceding discussion on species 
distribution, it is apparent that electrokinetic interaction 
between mercury ccsnplexes and the active carbon or coal surface 
is significant. Active carbon and coal surfaces are generally 
considered to bear net negative charges (Bean, Can^bell and 
Anspach, 1964? Pandey, 1978), Thus, for successful pr^iction 
of mercury c 3 oncentration-time profiles in a column system, 
the applied model should be incorporated with the contribution 
of electrokinetic interacticjn. However, this would make the 
model more <xmplex leading to difficulty in mathematical 
soluticxi. 
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Referring to the Famularo concept of micropore-fliacroshell 
particle (Famularo, Mueller and Pannu, 1980) it is recognized 
that the sizes of the different mercury conplexes are also 
in^ortant for soarption. The mercury complexes may fail to 
enter into the minute pores even if the potential for sorption 
and electrokinetic attractions are there. If the mercury 
species have restricted entry into the pores of the sorbent 
because of size consideration, the whole concept of sorption, 
electrokinetic interaction and diffusion is no more totally 
valid. 

B. Forms of Mercury — Inorganic, organic and Chelated Species 
The potential for microbial methylation of mercury by 
fungi and bacteria exists under aerobic as well as anaerobic 
conditions (Jemelov, Landner and Larsson, 1975), In a 
column containing the sorbent, the inorganic mercury ccxning 
in as feed might be transformed into organic mercury. This 
possibility brings in two problems, viz.^ a) the model 
selected for only inorganic mercury removal becomes insensitive 
and a new bacteriologically activated carbon model similar 
to that of Ying and Weber (1979) needs to be applied, and 
b) even though organic mercury removal is more than that of 
inorganic mercury by the same system, the advantage is only 
apparent because the small amount of organic mercury present 
in the effluent might be more hazardous compared to the 
higher concentrations of inorganic mercury in the column-feed, 
Predichion of the actual environment inside the column 
is rather difficult. Apparently, aerobic environm^t should 
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prevail in the coluirai for a continiK)us feed system as the 
influent has sufficient dissolved oxygen. On the other hand, 
active carbon through some mechanism has the ability to 
preferentially adsorb oxygen from air and oxygen deficioicies 
might occur even in the ccwnplete absence of biological acti- 
vity (Strudgeon and Lewis, 1980). If the dissolved oxygen 
content of the column influent and effluent are known, 
predictions abait the inside environment may be made. However, 
the possibility of existence of son^! suiaerobic regions cannot 
be ruled out even if the column effluent contains detectable 
levels of dissolved oxygen. 

The reported removal of Hg(ll) by either active carbon 
or coal system was not very high (Humenidc and Schnoor, 1974; 
Pandey, 1978). Thiem, Badorek and O’Connor (1976) treated 
aquecus solutions of ntercury with various combinations of 
active carbon, EDTA, tannic acid, citric acid and calcium. 

The results of their experiments indicated that the presence 
of chelating agents enhanced the removal of mercury fsrcxn 
water using povKiered active carbon as sorbent. Now if in 
column systems also chelating agents are added, the predicted 
concentration- time profile on the basis of carbon-only systen 
would be different. Moreover, it would be extremely diffi- 
cult ‘•o propose and/or apply a mathematical sorption model 
.incorporating the contribution of the chelating agents. 

C. Attainment of Equilibrium 

Extendol contacting periods should be allowed during 
isotherm evaluations to ensure that equilibrium is attained; 
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wherever possible, powdered active carbon isotherms should be 
used instead of granular carbon isotherms (Peel and Benedek, 
1980) . Isotherm data are very important for column-system 
effluent prediction. The available isotherm data for mercury 
(Humonick and Schnoor, 1974; Pandey, 1978) were not generated 
with the powdered form of the sorbent and it was not known 
whether true equilibrium had been attained because sanetimes 
it may take even five weeks to attain true equilibrium (Peel 
and Benedek, 1980). 

D. Analysis of Mercury 

Methods for analysis of total mercury have been a matter 
of toxicological concern and the subject of numerous reviews 
for many decades (Reimers, Burrow® and Krenkel, 1973). The 
methods have proliferated with the introduction of refined 
instrumentation, e.g,, atomic absorption, neutron activation. 
X-ray fluorescence and mass spectroscopy, and electrokinetic 
techniques. Low concentrations of trace metals can now be 
determined with precision. The different methods of mercury 
analysis reported in the literature ( Ballard and Thornton, 
1941; Fishman, 1970; Chau and Saitoh, 1970; Elly, 1973; 
Standard Methods, 1975; Jadcson and Dellar, 1979) are not 
discussed here, Greenberg, Thomas and Tamplin (1973) gave a 
table showing the comparison of some of the commonly used 
methods. The table, being considered to be an important 
reference for selection of ir^thods, is presented here. 
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Sanple 1 
( 0 pg/1) 

0,28 

0.50 

0.10 

1.00 

12.60 

Sample 2 
( 4 p.g/1) 

3.90 

4.60 

4.00 

1 .40 

9.60 

Sample 3 
( 6 \ig/l) 

5,40 

7.50 

5.80 

5.00 

9.00 

Sensitivity 

0.1 

0.3 

0,1 

0.5 

- 


It is apparent from the above brief discussion that the 

needed data on mercury sorption to be employed in mathematical 

models must be cnllected using the mercury analyzer. 

E. Hazards of Active Carbon and Coal as Sorbents for 
Mercury Removal 

It is interesting to note that a school of thought 
(Pendygraft, Schlegel and Huston, 1979) is against the wide- 
scale use of active carbon and charcoal as sorbents, even for 
organics removal for which these are so effective. They 


quoted Joshua Lederberg, a Noble Laureate, warning the U.S. 
National Drinking Water Advisory Council in 1976, "I have to 
tell you, you could make a horror story of the chemistry of 
what is in charcoal. If you were to take charcoal and 
disaggregate it and hydrogenate it, it would end up with a 
lot of horror substances that we would really seriously like 
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directed to the particular concern and I would most urgently 
and most seriously advocate them before v/e prcmilgate the 
use of charcoal adsorption on a large scale for our public 
water supplies." In an experiment, twentyfour metals were 
leached frcxn the Virgin granular active carbon, including lead, 
mer< 3 ury, cadmium and arsenic (Pendygraft, Schlegel and Huston, 
1979). Ihis leads to the belief that active carbon or coal 
is not without blemish as an adsorbent. 
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4. WEIOilNG MODIFIED /iPPROACHES t'fllTH iWAILf^LE DATA 

A review of the available sorption models on the basis 
of data need and preliminary assessment in terms of their 
applicability to mercury systems (Chapter 2) ccupled with a 
consideration regarding the special aspects related to mercuiy 
sorption led to the inclusion of the following approaches for 
final consideration. These approaches are based on the 
sorption models discussed in Chapter 2. 

Appiroach A 

Second order kinetics using Ihomas solution (Thomas# 

IS 44) and supplemented by Hiester and Vermeulen ( 1952) empl- 
oying (i) theoretical determination of the sorpticai rate 
constant according to Hiester and Vermeulen ( 1952) and ( ii) 
graphical determination of the sorption rate constant accor- 
ding to Miyamoto and Tram (1979). 

Approach B 

A modified approach in line with the lumped parameter 
method used by Keinath (1975) and coupling the diffusicai 
model of Newman ( 1931) , 

Approach A 

A notable limitation of the Thomas solution is that 
the rmdel is valid for Langmuir isotherm only. For Freund- 
lich and other isotherms# the model is not applicable, 
Miyamoto and Tram (1979) also developed their graphical 
method of predicrting the sorpticai rate constant for 
conditions satisfying Langrruir isotherm only. Also# 
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thGoretical detennination of the sorption rate constant 
according to Hiester and Vermeulen ( 1952) is accomplished 
using Langmuir isotherm data. The Thomas solution ( equati^i 
16) was applied to predict effluent concentration-time 
profiles in the mercury-coal syston of Pandey. ( 1978) and 
mercury-active carbon system of Humenicdc and Schnoor { 1974) . 
The Langmuir isotiierm plots for the two systems under consi- 
deration are shown in Fig. 3. 

Mercury-Coal system 

Equilibrium parameter, r, in equation (13) = ^ * 

o 

The calculated equilibrium parameter, r, is nearly unity for 
K = 0.3744 1/mg (Fig, 3b) and C = 10 jig/L (Pandey, 1978) 

W 

‘ k ^ ^ KSL 

Column capacity parameter, s s= (17) 


F or r " 1 , ^ — 




D. 


( 19a) 


the parameters were defined in equation (19). 'Ihe column 
capacity parameter was calculated frcsn the parameters shown 
in a tabular form in the following page. 

t ^ 

Ratio of capacity parameters, — = g — ^ — (14) 

G E 


For r 5= 1 and s = 4,75, t/s = 0,25 ( frcan Fig, 2) for break- 
throu'^h (C/C = 10 percent). To reach breakthrough, it was 

O ' 

calculated from equation (14), the time required would be 

■ 3 3 

about 27 months for = 1707.3 x lO and vf^ = 17.64 cm . 

The period of 27 months for reaching breakthrough is too 

long a time to accept eventhough the available data 
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P arameter Source of S valuation value 



cf. kpA^ confuted in 

Approach B 

0.032/hr 

°G 

*^co b 

C fr. 

O E 

1707.3 X 10^ 


Pandey ( 1978) 

0,078 on 

€ . 

Pandey (1978) 

0.5 


determined in Approach B 

-9 2 

2 X 10 cm /sec 


assumed 

0.5 

D 

pore 

equation ( 21) 

0,5 X 10“^ cm^/si 


equation ( 20) 

0.0109 549 /hr 

k 

equation ( 19a) 

0.03225/hr 

eb 

Pandey (1978) 

3 

0,7 g/cm 

L 

Pandey ( 1978) 

45 cm 


Pandey ( 1978) 

2 

8.15 ml/min/on 

K 

Fig. 3 

0.3744 1/mg 

•s 

Pig. 3 ■ 

12.195 mg/g 

S' 

equation ( 17) 

4.75 




c/s, mg/g C/S, mq/q 


1.4 r S 02.715172 rr.g/g 



C, mg/l 



FIG. 3 LANGMUIR ISOTHERM PLOTS FOR MERCURY 
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{Pandey^ 1978) shows that there was no sign of breakthrough 
up to run time of 8 days. 

Following the graphical approach of Miyamoto and Tram 
( 1979), the value of the rate constant was determined from 
Fig, 4(b) to be 0,068/hr, The calculated column capacity 
parameter, s, for this rate constant was 10,02, For r = 1 
and s s= 10,02, t/s = 0.45 ( f 3Com Fig. 2) for breakthrough 
(C/Cq = 10 percent). For breakthrcugh, the required calculated 
time (from equation 14) was 49 months which is, again, a 
totally unacceptable prediction, 

Mercury«-a.ctive Carbon System 

The equilibrium parameter, r, the column capacity para- 
meter, s, and the ratio of the capacity parameter, t/s, were 
similarly calculated for this systan also from the following 
available data (Humenick and Schnoor, 1974); K = 0,08855 1/mg 

(Fig, 3a) , S = 12.715172 mg/g (Fig. 3a), L = 107 cm, k * 

2 

0,05/hr (assumed), v^ =9,4 ml/min/cm , = 3 mgA, fg = 0.5, 

vf_ = 606,7 csn , The column capacity parameter, s, was found 

E , 

to be 2,42 and was calculated as 3840, the latter still 
being apparently veiy high. The capacity ratio, t/s, was 
found to be 0,1 for breakthrough (from Fig, 2). For attaining 
this breakthrough the required time, calailated from equation 
(14), was 17 days and 6 Ccanparing with the available 

data (himenick and Schnoor, 1974) it was seen that the 
predicted results were very encouraging as the breakthrough 
was not achieved in the actual column evoi after 5 days. 



2 3 

TIME (t), hr 

FIG. 4 PLOTS FOR DETERMINING THE RATE 
CONSTANT. 
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■The fact that the Ihomas solution was reasonably appli- 
cable to the me rcaary -active carbon system (C^ * 3 mg/1) and 
failed miserably in the mercury— coal system *= lOjig/1) 
demonstrates that the influent concentration is a very signi- 
ficant factor and that if it is very low, the approach becomes 
invalid. From the calculations it is also seen that a very 
low influent concentration gives an abnormally high value of 
Dq which ultimately increases the required run time for 
breakthrough. Eventhcugh that the break thrcsugh is not indep- 
endent of concentration is acceptable, the enoinnous increase 
in the theoretical run time required with decreasing concen- 
tration only demonstrates that the model is inapplicable in 
the case of mercury concentrations in the microgram range. 
Approach B 

A review of the Keinath Model (Keinath, 1975) shows three 
limitations, viz,, ( i) the model is not nondimensional, ( ii) 
the CSMP/360 program used for soluticai is not available with 
the dec-10 system at I.I.T,, Kanpur, and ( iii) this model is 
also limited to Langmuir isotherm conditicns only. 

Ihe first two limitations were eliminated by proposing 
a modification of the Keinath model in the nondimensional 
form and numerically solving the equations by the fourth 
order Runga-Kutta method. The following lumped x^rameter 
model equations developed by Keinath were discussed in 
Chapter 2. 
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dt " ^^n-1 ■ V 

( 11) 

^ - R 

dt ' A 

(12) 


Keinath ( 1975) discussed in details about the determin" 
ation of the rate of sorption, R^, and the rate equations for 
pore dif fusion,pore-surf ac3e diffusion and film diffusion were 
givon. Eventhough various theories including the penetration 
and free surface models have been advanced to describe trans- 
ient diffusion through an external boundary layer, the straight- 
forward relationship given in the following equation ia the 
only form that is generally applicable to simulation of the 
dynamics of a columnar adsorber 

= it = Vs - =*> 

where, k^, = film-diffusion controlled mass transfer 
coefficient, cm/hr, 

Ap = external interfadLal transfer area for 
sorbent, 

'k 

C s= solution-phase concentration of solute consi- 
dered to be in equilibrium with the outer 
surface of the sorbent particle, mole /L, 

“ t for Langnuir isotherm), and (30) 

%ax “ 

B = Langmuir energy term, lA^r. 
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Because of the difficulty in describing the hydrt^dynamic 
conditions in the boundary layer^ Keinath ( 1975) determined 
the mass transfer coefficient^ from a correlation of the 
mass transfer factor/ with Sherwood number, Reynolds 
number and Schmidt number as follows 

where, ss , mass transfer factor, 

^Sh “ Sherwood number, 

Ng^ = Schmidt number, 

, Reynolds number, 

U = U/A, velocity of flow, can/hr, and 
^ = kinematic viscosity of solution, an^/sec. 
The equation used by Keinath (1975) for determining the mass 
transfer factor, j^, was 

= 5.7(Nj^g)“^-'^® (32) 


The details of the formulation of the modified non- 
dimensional approach developed for application in the present 
study are given below 


Y 

n 


-S3- 

^ax 


(33) 


(34) 


T 



( 35) 
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E = 


C B 
o 


e C, 


G = 


^ ^raax 

evky 

u 


Fran equations (30) and (34) 

y 

b( 1 - Y) 


★ y 

c = 


From equations (33) , (36) and (39) 
E Y 


★ 

= 




( 36) 

( 37) 

( 38) 

( 39) 

(40) 


Equation (33) through equation (40) were substituted in 
equations (11) and (12) to obtain 


dZ 

n 

dT 


1 " 2„) - G(Z„ 
n— JL n n 


I 


E Y 



(41) 


E Y 


= GF( Z, 


n 


n 


n 


(42) 


The initial oDnditicMis: at T = 0, z = 1 and Y = O. ibe 
details regarding the different parameters used in this model 
are givai in Appendix B . 

The fourth order Runga-Kutta method was used for 
numerioal solution of the equations (41) and (42). The 
prograni listing for the solution is given in Appendix C. 

The coefficient of diffusivity of mercury conplexes in 
solution, is an important parameter v^ich is difficult to 
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determine (Section E; Chapter 2). Keinath ( 1975) did not 
explain how he calculated for the present study, was 
calculated on the basis of the diffusion model of Nevsnan 
( 1931), Frcan the batch kinetic data of Pandey (1978) and 
Fig, 1, was determined as 2 x 10 cm /sec for the mercury- 
coal system. The diffusion model-predicted concentrations 
for a batch sorption system on the basis of this value of the 
coefficient of diffusivity, and the experimentally deter- 
mined concentrations were plotted in Fig. 5, It was inter- 
esting to note that although the predicted concentrations were 
showing an initial higher rate of the decrease in the liquid- 
phase mercury concentration, after stxnetime ( 1-3 hr) the 
predicted concentrations were higher than the actual concent- 
rations, This showed that when the concentration reduces to 
a certain limiting value, in addition to diffusion there 
exists one or more controlling phenomenon. It could well be 
convection as in the case of the Convection Limited Model of 
Hendricks (1980) (Section K? Chapter 2). It might be 
diffusion-ccxivection instead of diffusicai-only , which causes 
greater ranoval of mercury from liquid phase than the predicted 
values. This itself shows that the coefficient of diffusivity, 
is difficult, if not almost impossible, to determine very 
accurately for the systems under consideration. Nevertheless, 
the value of was accepted (or rather asaimed) 

as 1 - 2 X 10*"^ csn^/sec for the mercury- coal system. 

8 2 

Similarly, 8 x lo“ on /sec was considered acceptable for 
the mercsiry -active carbon sys ton using the data of Ifiamenick 
and Schnoor ( 1974) . 
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FIG. 5 DIFFUSION MODEL PREDICTED AND EXPERIMENTALLY DETER- 
MINED CONCENTRATIONS IN MERCURY-COAL BATCH SORPTION 
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On the basis of the different parameters defined in 
Appendix B/ the numerical solution predicted effluent concen- 
tration-time profiles for each of the ten column elements are 
shown in Piof, 6, It is seen that for mercury-active carbon 
system (Pig, 6a)/the predicted breakthrough occurred after 
2 days and 7 hr whereas Humenicdc and Schnoor ( 1974) in their 
experiment could not get the breakthrough even after 5 days. 

For mercury-coal syston (Fig, 6b) ^ the predicted breakthrough 
occurred after 7 days whereas Pandey ( 1978) could not get the 
breakthrough even after 8 days. Unfortunately, data were not 
reported for run time beyond 5 and 8 days for the mercury- 
active carbon and mercury-coal systems, respectively and hence 
the degree of discrepancy could not be ascertained. However, 
it cxjuld be very well said that the modified lumped parameter 
approach predicted an early breakthrough, i 

The fact that the approach considered included several 

significant parameters must be kept in mind while any attempt I 

i 

is made to criticd.ze its apparent inadequacy per se in terms 
of predicting the performance of a real-world system. Para- i 

' • ' • - ■ i 

meters like the rate coefficient of sorption, diffusivity, | 

i 

coefficLait of mass transfer, etc. , are difficult to determine | 

j 

without reservations. Now, if an approach fails, whether | 

the blame goes to the conceptual formulation of the approach 
or the limitations involved in determination of individual 
parcimeters, is a matter for open discussion. The approach 
v^as based on the lump^ parameter method vhich is widely 
used in process dynamics (Priedly , 1972) and has been 
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FOR MERCURY -ACTIVE CARBON AND 
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applied by other workers to the sorbent column-system (Wilson 
and Clarke/ 1979), So one has to make sure that the signi- 
ficant parameters are evaluated in a reliable manner before 
blaming the approach. 

Lack of reliability in evaluating the mass transfer 
coefficient/ kp/ and the effective area of transfer/ h, 
prevents precise weighing of the Approada B. One way of 
making an acid test of the successful applicability^' of the 
approach would be as follows; ( i) for a particular system like 
the mercury “Coal or mercury— active carbon system, experi- 
mental data are to be collected for different sets of influent 
concentration, velocity of flow, etc,,(ii) for one such set, 
the value of kpA may be determined using the model and other 
known parameters required, ( iii) this k^A may be considered 
to be constant for the particular system, ( iv) with this 
value of kpA and other known parameters, the concentration- 
time profile may be predicted for other sets of influent 
concentration, velocity of flow, etc., and (v) comparison of 
this predicted profile with the experimental one would give 
a good idea about the success of the approach, 

Keinath ( 1975) reported the predictai effluent concen- 
tration-time profiles only and did not include a ccanpariscsa 
with the experimental profiles. Also, he did not indicate 
the CPU time required for running the program in the computer. 
Thus v^.ether Approach B could give better results than the 
original lumped parameter method of Keinath (1975) solved 
by CSMP/360, could not be ascertained. Eventhough the 
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concept of the lumped parameter model as such might be 
without blemish, the fact that the model contains so many 
parameters is its main drawback. 

It is important to recognize that the set of differential 
equations (Equations (41) and (42)) are extremely stiff, i.e., 
the time constant for the liquid-phase equation is small in 
contrast to the time constant for the solid-phase equation. 
Numerical solution the Runga-Kutta method required 
specification of an integration interval that provided for 
the stable and accurate integration of the liquid-phase mass 
balance equation which responded most rapidly. As during 
the examputation it was seen that the increase of the integra- 
tion interval beyond 0,10 and 0,004. for mercury -coal and 
mercury-active carbon, respectively, gave rise to instability, 
the computations were done with these intervals which 
ultima tX3ly required pretty high CPU times. 

In spite of all these minus points, it should be 
recognized that the Keinath model is very useful as a starting 
point for any modified approach. Keinath (1975) himself 
hoped that the model would provide for the development of a 
generalized descriptive model which can be employed as a 
toolfotr understanding the ccxnplicated phencamenon of sorption. 
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5, SUMMARY AND CONCLUSIONS 

A good number of sorption models are available in the 
literature, most of them being related to organics sorption. 
None of these models were proposed for or applied to heavy 
metal sorption. To facilitate easy predictions of the removal 
of mercury by active carbon or coal in a column system. 
Investigations were made to find a reliable and handy choice 
from the listing of the available models. The available 
information on mercury species distribution in solution was 
studied and the importance of electrokinetic interaction, 
that had so far been neglected, was recognized. It was 
thought that the sizes of different mercury complexes should 
also be investigated into for predicting to vhat extent these 
woild enter the micropores of the sorbent, Ihe preliminary 
screening, based mainly on the availability of suitable data, 
eliminated all but only a couple of the available models, 
viz., the Thomas solution and the Keinath model. The latter 
was modified to the nondimensional form and as the library- 
program originally used was not available, the equations were 
solved by the fourth order Runga-Kutta method. Also, an 
approach based on the second order kinetics using Thomas 
solution and employing a theoretical as well as a graphical 
determination of the sorption rate constant was included. 
These constituted the major porticai of the investigation. 

It was recognized that the crux of the prcblon, which did 
not permit satisfactory testing of the approaches considered. 
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was inability to reliably evaluate the significant parameters 
involved . 

Considering the fact that the present study is a step 
forward in the area of modeling heavy metal sorpticai/ it is 
not possible to totally accept any available model or propose 
a now model* It may be concluded without hesitation, however, 
that none of the available organics sorption, models is 
acceptable for predicting mercury removal by active carbon or 
coal without major modifications. 
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6* SUGGESTIONS FOR FURTHER W5RK 

1, The important parameters like electrokinetic inter- 
acticai as well as species distribution an3 the size effects 
of the different probable mercuiiy cc3iT5)lexes in solution shaild 
be incorporated into any model for better simulation of mercury 
removal, in spite of the coitputaticnal ccxnplicacy that would 
be brought in "by such incorporation. 

2, Any available model limited to Langmuir isotherm 
conditions only should be modified for applicatiai to other 
common isotherm conditions like Freundlich, BET, etc. 

3, Appreciating that the sorption rate constant for 
Thomas solution of the second-order kinetics is of utmost 
importance^ efforts should be made to evaluate the rate cons- 
tant on the basis of more reliable experimental data, 

4, The liquid-film mass transfer coefficient may be 
calculated by the application of the modified Keinath model 
(Approach B) with the available sorption data on the basis of 
trial and error tuning the conputer. Once the transfer 
coefficient is found for a particular system, it may be used 
to predict effluent ccncentratiort-time profiles for the same 
system , 
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APPENDIX A 

LIST OF MATHEMATICAL SYMBOLS 


The following symbols were used 


A = 
P 


a = 


ap = 


^ = 



c, .c . 

So’ 


D 


D = 


= 


D, = 


'1 

pore 
£ 

j^E 


■I 
H( 2 ) 

I 


f. = 


external interfacial transfer area for sorbent, 
cm /g, 

i 

original amount of sorbate, mg/1, 

effective area for mass transfer for case of 

2 

fluid“phase transfer controlling, cm , I 

! 

effective area for mass transfer for case of | 

2 i 

transfer inside solid particles, an , 

Langmuir energy term, 1/hr, 

i 

influent concentration of solution-phase in a 
column system or initial concentration of I 

i 

solution-phase in a batch system, mole/1 or mg/l, 

solution-phase concentration of solute, rag/1, I 

j 

solution-phase ccncentration of solute at i 

equilibrium, mg/1, 

2 , 

molecular diffusivity, cm /sec, 

limiting distribution coefficient, | 

* 2 j i 

diffusivity of solute in solvent, cm /sec, i 

sorbent particle diameter, cm, | 

2 ' 

pore diffusivity, cm /sec, | 

■ . ■ I 

external porosity, 
internal porosity, 
error function, 

modified Bessel function of the first kind, 

! 

■ . ' ■ ' i 



57 


= mass transfer factor, 

k S5 sorption rate constant, hr"”^, 

kp = mass transfer coefficient for case of fluid- 
phase transfer controlling, cm/hr, 

kp = mass transfer coefficient for case of transfer 
inside solid particles, c3n/lir, 

K = Langmuir isotherm constant, 1/mg, 

“1 

K‘ = sorption rate constant, hr , 

L = length of column, cm, I 

N = sorption capacity coefficient, | 

N„ = Reynolds number, i 

N„ = Schmidt number. 

Sc I 

N„. = Sherwood number, i 

Sh I 

n = an integer, ! 

n = number of elements in the lumped parameter model, ! 

! 

Q = ultimate sorbent capacity, mg/g, 

« . 

q = solid-phase concentration of solute, mg/g, 

a s= solid-phase concentration of solute at equili- 

CO " 

brium, mg/g, 

r = sorbent particle radius, an, 
r = equilibrium parameter in 'Thomas solution, 

R = volumetric flow rate of fluid through fixed 

sorbent, ml/min, J 

r= rata of sorption, mole/g-hr, 

S = ultimate capacity of dry sorbent, mg/g, 

s = amount of sorbate sorbed per uftit weight of 

sorbent at saturation, mg/g. 
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s = column- capacdty parameter in Thomas solution, 
t = time, hr, 

t = solution capacity parameter in Thomas solution, 
U = solution volumetric flow rate, 1/hr, 

U £= velocity of flow, cm/hr, 

V = total volume of feed in column in Thomas 

solution, ml, 

V = volume of each element in lumped parameter 

model, 1, 

V = batch volume of liquid in Famularo model, 1, 

V = approach velocity, l/hr, 

V = bulk volume of column under consideration, ml, 

2 

= macroscopic water flow velocity, ml/on /hr, 

X 3 ? = amount of sorbate sorbed per unit weight of 

^ t 

sorbent, mg/g, 

x^ = bed depth in Logistic curve model, on, 

V X = amount of sorbate sorbed per unit weight of 
^cx>" eqm . 

sorbent at saturatim, mg/g, 

X , X 2 = amount of sorbate sorbed per unit weight of 

sorbent at time t^^ and t 2 respectively, mg/g, 

z = axial distance, cm, 

= density of dry sorbent in column, g/on , 
e = solid-solution ratio, mg/1, and 

^ = kinematic viscosity of solution, csn /sec. 



59 


CARBON 

A 

EP 

RHO 

EPSI 

QMAX 

B 

CO 

DL 

AREA 

U 

E 

P 


APPENDIX B 

PARAMETER INPUTS: MODIFIED KEINATH MODEL 



Coal 

Active- Carbon 


(Pandey, 

1978) 

(Itimenick and 
Sciinoor, 1974) 

mass o£ carbon/coal in 
bed, g 

24.74 

550,0 

external transfer area 
for cjarbon/coal , on^/g. 

35.00 

150.0 

diameter of carbon/coal 
particles, cm 

0.078 

0.118 

packed bed density, g/1 

700,0 

453.23 

packed bed porosity 

0.50 

0,45 

Langmuir ultimate uptake 
capacity, moles/g 

6.098X10”^ 

6.357x10“^ 

Langmuir energy term, 
1/mole 

0.749x10^ 

0.177x10^ 

solution-phase conccait- 
ration of solute, mole/1 

O.SXIO"*^ 

1.5x10’*^ 

diffusivity of solute, 
omVsec, 

1.22x10'^ 

-8 

8.0x10 

cross-sectional area of 
column, cm2 

0.785 

11.341 

solution volumetric flow 
rate, 1/hr 

0,384 

6.396 

nondimensicaial parameter 

267.094 

3.764 

nondimen sional parameter 

0.586xl0”® 

0.234x10’’^ 

nondimensional parameter 

8.422 

631.285 


G 
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APPENDIX C 


COMPUTER PROGRAM 


C MODIFIED KEINAIH MODEL 


C 


C 

C 

C 


8 

C 


1 


2 

G 

100 


; 10 




Meraary-Coal System 

REAL K0,K1,K2,K3 

DIMENSION X{ 20) ,Xl( 20) ,X2( 20) ,X3( 20) ,FX( 20) , 

1K0( 20) ,K1( 20) ,K2( 20) ^K3( 20) 

COMMON E,F^G 

OPEN( UNIT=1 ,FILE= • MPP2 .RES ’ ) 

DATA E/.2670940E+03/,F/.5857157E-06/,G/.84216937E+01/ 

DELT=0,10 

KKFsSOOO 


INITIALIZATION 


T=0,0 

DO 8 1=1/20 
X( I) =0.0 
KK=0 


TYPE 1 
WRITE(1/1) 

F0RMAT(49X/ 'KEINATH MODEL MERCURY»C0AL SYSTEM’ //42X/ 

1' CONCENTRATION PROFILES IN SIMPLE PACKED BED COLUMN’/ 
242X/50( *=’ )) 

TYPE 2,E/F/G/KKP/DELT 
WRITE( 1/2) /E/F/G/KKF/DELT 

FORMAT (//1X/‘E=’ /E14 .7 /2X/ ’ F=’ /El4 .7,2X, ’G=’ /E14.7/2X, 
l'KKF=' /I6,2X/’DELT=' ,F5.3) 


CALL DERIV(X,FX) 

DO 10 1=1/20 
K0(I)=FX( I)*DELT 
Xl( I)=X( I)+0.5*K0( I) 

CALL DERIV( X1/FX) 

DO 11 1=1/20 
K1(I)=FX(I)*DELT 
X2( 1) =X( I) +0.5*K1( I) 

CALL DERIV(X2/FX) 

DO 12 1=1,20 
K2(I)=FXCI)*DELT 
X3(I)=X(I)+K2( I) 
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CALL DERIV{X3,FX) 

DO 19 1=1,20 

19 K 3( I ) =FX( I ) *DELT 

DO 25 1=1,20 

25 XC l)=X(l)+(K0(l)+2.*Kl(l)+2.*K2(l)-f^C3(l))/6.0 

T=T+DELT 
KK=jac+l 
SUMQ=Q,0 
DO 60 1=11,20 
SUMQ=SUMQ+X( I) 

60 CONTINUE 

if(kk.ne.ki<;f) go to loo 
TYPE 5,T,SUMQ,(X( I) ,1=1,20) 

WRITE( 1,5) ,T,SUMQ,(X( I) ,1=1,20) 

5 FORMAT(/lX, 'T=‘ ,Fll .3, 2X, ' SUM Y( N) =* ,E15 .7/lX, ' Z’ ,10El3.7y 

llX,*y* ,10E13.7) 

KK=0 
PAUSE 
GO TO 100 
STOP 
END 

SUBROUTINE DERIV(U,DU) 

DIMENSION U( 20) ,DU( 20) ,PC 10) 

COMMON E,F,G 
DO 10 1=1,10 

10 P( I)=G*(U( I)-(E*U( I+10)/( 1.0-U(I+10)))) 

DU(1 ) =( 1 .0-U( 1 ) )-P( 1 ) 

DU( 11)=F*P( 1) 

DO 20 1=2,10 
DU(I)=(U(I-1)-U(I))»PC I) 

DU( I+10)=F*P( I) 

20 CONTINUE 
RETURN 
END 


C 



